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Abstract 

For random percolation at pc, the probability distribution P{n) of the number 
of spanning clusters (n) has been studied in large scale simulations. The 
results are compatible with P{n) ~ exp(— an^) for all dimensions. We also 
study the variation of the average size (mass) of the spanning clusters when 
there are more than one spanning cluster. While the average size of the 
spanning clusters scales as usual with where D = d — (3/u for any number 
of clusters, it shows a smooth decrease as the number of spanning clusters 
increases. 
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More than one spanning cluster (when the spanning is considered in one direction) at the 
percolation threshold for dimensions below 6 have been shown to exist in a number of recent 
studies [1-6]. The numerical evidence for nonunique spanning cluster was already there in 
five dimensions in de Arcangelis' study of spanning clusters [1] nearly a decade ago, which 
was, however, interpreted as a finite size effect. 

Recently, there have been some controversies regarding the behaviour of the probability 
distribution P{n) of the number of spanning clusters (n). Using mathematical arguments, 
Aizenman has shown [2,5] that P{n) ~ exp(— an^) in two dimensions. Aizenman [2] also 
conjectures that in any dimension P{n) ~ exp(— an^). In simulations reported earlier 
[3], we had remarked that the distribution appears to be a simple exponential (for 4 and 5 
dimensions, only for which considerable data was available). But for reasons given below 
and on the basis of our present results, this now seems less clear. In the earlier study ([3], 
referred to as I henceforth) , we obtained these probabilities by varying the site concentration 
p (the initial value taken close to the known value of Pc) until a spanning cluster occurs and 
then counted the number n of such clusters. Here we use the p^ for the infinite lattice thereby 
increasing the efficiency of the program. However, the numerical results arc then different. 
We also later realised that the averaging in the previous data in 1 was actually being taken 
over a number of configurations lesser by a factor of 136 than the intended number due 
to an error in the parallelisation of the program. However, the main results of I have been 
verified to be qualitatively correct; the probability of getting more than one spanning cluster 
does remain finite. In fact, the present definition indicates a higher probability with lesser 
error and as an example we have presented the result for two dimensions in fig. 1 for this 
probability. Here, as in I, the probabilities have been normahsed by a factor of (1 — P(0)), 
which is the spanning probability. 

Here also, we have simulated hypercubic lattices with L'^ sites in d dimensions with 
helical boundary condition. The number of configurations over which averaging is done 
varies from 10"^ to 10^ according to the sizes of the lattice. Simulations are carried on for 2 
to 5 dimensions. 

We have plotted (fig. 2 and 3) the logarithm of the probabihties against (a) n and (b) 
ri^ for dimensions 2 to 5 and also against (c) n'^ for dimensions 4 and 5. It is observed 
that plot (b) gives a straight line over a wider range of n in comparison to (a) and (c), the 
latter correspond to a simple exponential and the exponent conjectured in [2] respectively. 
We also get support for an exponent equal to 2 by a different analysis. Assuming that 
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P{n) ~ exp(— an'''), let 



z = ln(P(n))/ ln(P(n + 1) ~ {n/{n + 1))^ (1) 

We have first plotted double-logarithmically z against to see whether we really get 
a straight line. This has been done for five dimensions only, where we have obtained the 
maximum number of spanning clusters. It is difficult to determine from fig. 4a which straight 
line of the three (with slopes = 1.00, 1.25 and 2.00) is the best fit to the data (for 1). 
This is compatible with the fact that there are some regions in each of the figures 2a, 2b, 
3a, 3b and 3c where a straight line can be fitted. However, it is also true that the data 
represented by the □ and x in fig. 4a, have a slope close to 7 = 2, for which the we have 
the best statistics (L = 15 and L = 13). 

We have also calculated x = ln(z)/(ln(n/n + 1)) for different sizes which should approach 
7 forn — > cxD according to (1). In two dimensions, we do not have enough clusters to calculate 
this quantity. In order to check whether the values are affected by the presence of a prefactor 
(which has been neglected in (1)), we have also calculated x with the absolute values P{n) in 
eqn (1) replaced by the ratios P(n)/(1 — -P(O)). We find that x approaches 2 in dimensions 
3 to 5. The results are shown in fig. 4b-d, where we show the values of x for both cases. 
The o represents the case when the prefactor is included. For three dimensions, in fig. 4b, 
the data close to x — 6 correspond to n = 1 and the one close to 2 to n = 2. With no 
prefactor, the data is above 2 for n =1 and slightly below 2 for n = 2. For d = A and 5 (fig. 
4c and fig. 4d), the values of x approach 2 from below as n is increased when there is no 
prefactor and from above when there is a prefactor (in general). We notice that the values 
of X from the two different measures come closer as n is increased in all dimensions and it 
also approaches 2 clearly in 4 and 5 dimensions where we have data for a larger number of 
clusters. In fig. 4c, data for n — 1,2 and 3 are shown while in fig. 4d, data for n — 2,3,4,5 
and 6 are shown. Although there are some fiuctuations in the data, apparently the exponent 
is closer to 2 than to 1. Our conclusion is therefore that for intermediate n the probabihty 
distribution behaves as a stretched exponential P{n) — exp{—a-nP'), with the exponent close 
to 2 in all dimensions. Apparently, the value of a decreases with dimensionality (~ 0.5 for 
(i = 4 and ~ 0.2 for d = 5). 

We have also studied the sizes or masses (i.e., the number of occupied sites in the cluster) 
of the spanning clusters. The average size of the spanning clusters indicates that even when 
there are more than one cluster, the size still scales as L^, with D — d — P/u, as in the 
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case of unique cluster [7] . However, there is a decrease in the average size as more spanning 
clusters appear and this is also expected. In two or three dimensions, where we have at the 
most 2 or 3 clusters respectively in general, the average size of spanning clusters versus the 
sytem sizes L data is shown in fig. 5a. For 4 and 5 dimensions, where more clusters are 
obtained, the scaled data (m/L^), where m is the average mass, against the cluster number 
is shown in a log-log plot (fig.Sb). It is difficult to verify whether the variation with n is a 
power law because of the fluctuations for large number of clusters. These fluctuations are 
unavoidable as the averaging for cluster sizes for larger n is done over a smaller number of 
cases (as the probabilty decreases with number of clusters). However, the smooth decrease 
is clearly indicated. 

We also calculated the second and third moments of the masses where the qth moment 
is defined as = ^ ~ "^Y^ to get an idea of the probabihty distribution of the 

sizes. While mij-fr? remains more or less a constant for different system sizes in five dimen- 
sions for six clusters, m-ijfr? remains positive and fairly constant indicating an asymmetric 
distribution : possibly one large cluster and several smaller clusters exist. 
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Stauffer for suggestions, discussions and encouragement. She also thanks A. Aharony and J. 
P. Hovi for discussions. Acknowledgements are also due to German Israeh foundation and 
to HLRZ Julich for 10^ hrs on Paragon. 
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Figure Captions 



Fig. 1 The variation of the probabihty of getting more than one spanning cluster is shown 
against different system sizes in two dimensions. 

Fig. 2 Probabihty distribution {P{n)) of the number of spanning clusters (n) for 2 {o L — 
150) and 3 (upper data points) dimensions against (a) n and (b) n^. The different symbols 
are for sizes : + L = 100, □ L = 125, x L = 150, A L = 175, * L = 200 in 3 
dimensions. 

Fig. 3 Probabihty distribution (-P(n)) of the number of spanning clusters (n) for 4 (lower 
data points) and 5 (upper data points) dimensions against (a) n and (b) ri^ and (c) against 

d 

nd-i. The different symbols correspond to different system sizes : o for L = 21, + for L = 
19, □ for L = 17, X for L = 15 and A for L = 13 for d = 5 and * for L = 39, filled □ for L 
= 35, filled A for L = 30, o for L = 27, and + for L = 23 for d = A. 

Fig. 4a The values of z against as found for 5 dimensions. The different symbols 
represent different sizes : oL = 21, + L = 19, □ L=15 and x L = 13. The three 
straighlines have slope = 1.00, 1.25 and 2.00 (from left to right). 

Fig. 4b The values of x vs. L as found for 3 dimensions. The o arc for the ratios and the 
+ for the absolute values of P{n). n increases as the value of x decreases. 

Fig. 4c The values of x vs. L as found for 4 dimensions. The o are for the ratios and the + 
for the absolute values of P{n). n increases as the value of x decreases for o and vice versa 
for the +. 

Fig. 4d The values of x vs. n as found for 5 dimensions. The o are for the ratios and the + 
for the absolute values of P{n). The data correspond to different sizes L=13toL = 21. 

Fig 5 a The variation of the average mass of the spanning clusters against the system sizes 
ioT d = 2 (o for no. of spanning clusters n = 1 and + for n = 2) and d = 3 {O ior n — 1, x 
for n = 2 and A for n = 3). The system sizes in d = 3 have been multiplied by 10 in order 
to be shown in the same range. 

Fig 5b The variation of the scaled average mass (fh/L^) of the spanning clusters against 
the the number of clusters in 4 and 5 dimensions. The average masses in four dimensions 
have been divided by 2 for better viewing. The different symbols correspond to different 
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system sizes : o for L = 21, + for L = 19, □ for L = 17, x for L = 15 and A for L = 
13 for d = 5 and * for L = 39, filled □ for L = 35, filled A for L = 30, o for L = 27, 
and + for L = 23 for d = 4. 
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Scaled average mass of spanning clusters 
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